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Abstract: In the paper, we consider a multi-dimensional bipolar hydrodynamic model 
from semiconductor devices and plasmas. This system takes the form of Euler-Poisson 
with electric held and frictional damping added to the momentum equations. We show 
the global existence and LP convergence rates of planar diffusion waves for multi-dimensional 
bipolar Euler-Poisson systems when the initial data are near the planar diffusive waves. A 
frequency decomposition and approximate Green function based on delicate energy method 
are used to get the optimal decay rates of the planar diffusion waves. To our knowledge, the 
LP{p G [2,-|-oo])-convergence rate of planar waves improves the previous results about the 
L^-convergence rates. 
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1 Introduction. 

In this paper, we consider the following bipolar Euler-poisson system (hydrodynamic 
model) in three space dimension: 


dtp~^ + div(/9’''ri''') = 0, 

dtip'^uf) + diY{p+ufu+) + dxiP{p'^) = -p'^uf -h p'^d^iCj), 1 < i < 3, 


< dtp~ + div(/9“ti“) = 0, 

dt{p~u~) + div{p~u~u~) dxiP{p~) = -p~u~ - p~dxi(t), 1 < i < 3, 
. Ac/. = p+ - p-, lim|3,|^^ |V(/>| = 0, 


( 1 . 1 ) 


with initial data 


(p=^,w=^)(x,0) = {p^{x),u^{x)), 


( 1 . 2 ) 
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where are the two particles’s densities, = {p^u^, p^u^, p^u^) are current densities, 

(j) is the electrostatic potential, and P{p^) are pressures. As usual, we assume the pressure 
P{p) be smooth function in a neighborhood of a constant state p* with P'{p) > 0. The bipolar 
Euler-Poisson equations are generally used in the description of charged particle fluids, for 
example, electrons and holes in semiconductor devices, positively and negatively charged ions 
in a plasma. This model takes an important role in the fields of applied and computational 
mathematics, and we can see more details in na 1211125] etc.. 

Due to their physical importance, mathematical complexity and wide range of applica¬ 
tions, many efforts were made for the multi-dimensional bipolar hydrodynamic equations from 
semiconductors or plasmas. Li |18] showed existence and some limit analysis of stationary 
solutions for the multi-dimensional bipolar Euler-Poisson system. Ali and Jiingel [I], Li and 
Zhang [16] and Peng and Xu [23] studied the global smooth solutions of the Cauchy problem 
for multidimensional bipolar hydrodynamic models in the Sobolev space P[^(K‘^){1 > 1 + ^) 
and in the Besov space, respectively. Ju m discussed the global existence of smooth solutions 
to the initial boundary value problem for the three-dimensional bipolar Euler-Poisson system. 
Li and Yang m and Wu and Wang m showed global existence and decay rate of the 
smooth solutions to the three dimensional bipolar Euler-Poisson systems when the initial 
data are small perturbation of the constant stationary solution. Huang, Mei and Wang [7] 
showed large time behavior of solution to n-dimensional bipolar hydrodynamic model for 
semiconductors when the initial data are near to the planar diffusion waves. Ali and Chen 
[2] studied the zero-electron-mass limit in the Euler-Poisson system for both well- and ill- 
prepared initial data. Lattanzio m and Li m investigated the relaxation limit of the 
multi-dimensional bipolar isentropic Euler-Poisson model for semiconductors, respectively. 
Ju, etc. [13] discussed the quasi-neutral limit of the two-fluid multi-dimensional Euler- 
Poisson system. Moreover, it is worth to mentioning that there are a lot of reference about 
the one-dimensional bipolar Euler-Poisson equation, and the interesting reader can refer to 
[3[ (5] m m m m dni mi mi mi mi and the reference therein. In particular, motivation by 
Gasser, Hsiao and Li [4] found that the frictional damping is the key to the nonlinear 
diffusive phenomena of hyperbolic waves, and investigated the diffusion wave phenomena of 
smooth “small” solutions for the one-dimensional bipolar hydrodynamic model. Huang and 
Li [B] also studied the large-time behavior and quasi-neutral limit of L°° solution of the one¬ 
dimensional Euler-Poisson equations for large initial data with vacuum. That is, they showed 
that the weak entropy solution of the one-dimensional bipolar Euler-Poisson system converges 
to the nonlinear diffusion waves. Then Huang, Mei and Wang [7] showed the planar diffusive 
wave stability to n{n > 2)-dimensional bipolar hydrodynamic model for semiconductors, and 
obtained the optimal and L°° decay rates. In this paper, we are going to reconsider global 
existence of the smooth solution for the multi-dimensional bipolar Euler-Poisson systems, in 
particular, we try to establish the U’{p G [2, -|-oo]) convergence rates of planar waves. 

In the following discussion, we assume that the initial data are a small perturbation of 
the diffusion profile constructed later with small wave strength. Let the initial data p^{x) be 
strictly positive and satisfy 

liin = p±, 

where p± > 0 are two far field constants with p- ^ p^. Similar as the consideration of planar 
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diffusion waves of damped Euler equations in |28t [20] , to define the multi-dimensional planar 
diffusion wave, we first consider the one dimensional diffusion equation 

dtw = P{w)xixx, (1-3) 

which can be derived from the bipolar Euler-Poisson equations with the relaxation terms in 
one dimensional case by imposing the Darcy’s law, cf. Then a multi-dimensional 

diffusion wave w{x,t) is a one dimensional profile in multi-dimensional space. That is, 
w{x,t) = W (xi/Vl + 1) is a self-similar solution of the equation (11.311 connecting two end 
states p± at xi = Too. Denote C = , then W{Q) satisfies 

= d(^{P'iW{C))dcW). 

Eor simplicity, let the initial velocity u^{x) vanish as xi —>■ Too, that is, 

lim w^(x) = 0, 

xi^ioo ^ ^ ^ ’ 

which implies that there is no mass flux coming in from xi = Too. This assumption could 
be removed in a technical way similar to the argument for one dimensional problem because 
the momentum at xi = Too decays exponentially induced by the linear relaxation terms. 

We now recall the bipolar Euler-Poisson systems O) in one space dimension: 

dtp'^ + dx:^{p-^uf) = 0, 

dt{p'^uf) + dxiip'^ufuf) + dxiP{p'^) = -p'^uf + p'^E, 

< dtp~ + 9x1 {P~'H ) = 0> (1-4) 

dt{p~u{) + dxi{p~u{u{) + dxiP{p~) = -p~u{ - p~E, 

, dxiE = p+- p~,lim^j^^_ooE{xi,t) =0. 

Denote the solution of (II.4p by {p^^,uf, E){xi,t). When 

lim p^{xi,0) = p±, liin uf(xi,0)=0, 

Xl^±OD Xl^±OD 

the time-asymptotic behavior of {p'^,uf, p~,u]^){xi,t) has been studied in [Ij, which is shown 
to be a nonlinear diffusion profile governed by Darcy’s law. Roughly speaking, the solution 
p^{xi,t) converge to a same diffusion wave W (xi/Vl + 1) up to a constant shift in xi. Note 
that more detailed assumptions on the initial data of the one-dimensional problem (jl.4p will 
be specified in Theorem 2.1. 

In this paper, we will generalize this time asymptotic behavior towards a planar diffusion 
wave to three-dimensional case and establish the related {2 <p < oo) convergence rates. 

As in the consideration of planar diffusion waves of damped Euler equation in [281 [20] 
and of the bipolar Euler-Poisson system in [7|, we do not directly compare the solution of 
the problem (jl.ljl with the diffusion wave W{xi/y/l -|- t), instead, we will compare it with 
the solution of one dimensional problem (II.4j) . Eor this, without loss of generality, let us first 
assume the initial density p^(xi,0) in (II.4p satisfy 

r+oo 

/ (p=^(xi,0) - TT(xi))(ixi = 0. (1.5) 
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For the multi-dimensional problem, the shift function So(x') , where we used the notation 
x' = (x 2 ,X 3 ), can be chosen as in [28l[20] such that the initial density function satisfies 


r*+oo 


(/9^(x,0) -W{xi + (x')))dxi = 0 


Note that 5^(x') is then uniquely determined by 

(5j(x') =- / (p=^(x,0) - lF(xi))dxi, 

p+ - p- J-OO 

for 7 ^ /?_|_. Moreover, we assume that basically the shift is uniform in directions other than 
xi at infinity, that is, 

1 ^-1-00 

lim - / {p'^{x,0)— W{xi))dxi = df, 

Ix'l^-l-oo /9+ — P- 


Note that this assumption simplifies the problem and it remains unsolved for general pertur¬ 
bation when this assumption fails. An immediate consequence of this assumption is that 


lim 

|a;'|—>+oo 


^±(x') = 5: 


± 


And for simplicity, we assume <5^ = <5* be same constants. With these notations, the main 
purpose here is to show that the solutions (/9^,n^) of (jl.ll) converge to {p^,u^) with certain 
time decay rates, where 

p^{x,t) = p^{xi + 6{x' 

u^{x, t) = (nf (xi -I- 6{x', t),t), 0,0), 

_ _ ( 1 - 6 ) 
E{x, t) = {E{xi + 5{x', 0), 

5{x',t) = (5* -I- e“*((io(x') - 5*), 

in which p^,uf^ and E are solution of (jl.4p . In the following discussion, we will also assume 
that the shift generated by the initial data satisfies 


\d^A4ix')-^*)\<cii + \xf)-^, 


(1.7) 


for any multi-index /3 and any positive integers N. Here, C is a constant depending only on 
j3. This assumption implies that the shift So(x') decays to <5* almost exponentially. Again, 
this assumption can be reduced to the constraint on the initial perturbation. More precise 
construction of the background planar diffusion wave (p^,'U^)(x, t) and its properties will be 
given in Lemma 2.3 below. 

Throughout this paper, we denote any generic constant by C. The usual Sobolev space 
is denoted by IF®’P(R"), s G Z+, p G [l,oo] with the norm 

S 

||/||n.„ := ^ iiavilB, 

| a |=0 
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where 5“ used for 5“ without confusion. In particular, = R^(R"'). Set 


V^{x,t) = p^{x,t) - p^{x,t), 

U^{x,t) = {uf{x,t) - uf{x,t),uf{x,t),uf{x,t)), 

S7If = SJ(j) — E ( Note div(V0 — E) =0), 

K{x,t) = V~^{x,t) — V~{x,t), 

and also denote 

/ Xl 

V^{xi,x' ,Q)dxi, iyi^{x,0) = / x', 0)dxi. 

-OO j —CO 


( 1 . 8 ) 


(1.9) 


Note that the time derivative on the initial data can be defined by the compatibility of the 
initial data through the equation (| 2 .ip . 

Now we state the main results in this paper. Note that we consider only the spatial 
dimension n = 3 in this paper. However, we will still use notation n in the below theorem 
for convenience to extend our result to other high dimensional cases, since other higher 
dimensional cases can be considered similarly. 


Theorem 1.1 Let {p^,u^){x,t) in 11.61) be planar diffusion waves with a shift d{x',t) con¬ 
structed above. (See more precisely its properties in Lemma 2.3.) For k > A, assume that the 
initial data (p^,u^)(x, 0 ) satisfy the smallness assumption 

\p+- P-\ + + II(p=^ - P=^)(-,0)11^1 + \\{p^ - p^,u^ - u^)f,0)\\Hk < eo, 


where cq > 0 is a sufficiently small constant. Then 

(i) (Global existence) There exist unique global classical solution (p^,V(/)) to the system 
(Ufi-mM) that 


V^{t,x),U^{t,x) G C7([0,oo),R*=(R”)) n 67^(0, cx)), R*=-^(R")), Vf G 1E*^’®(R"). 
(a) (LT convergence) Moreover, for I 7 I < A: — 2, p G [2, 00 ], we have 

||a 2 H±|Up< Ceo(l + t)" 2 ('-F)-—, 

-.,1 n /I 1 \ I 7 I +2 

\\d2U^Lp< Ceo{l + t)-^^^-p^-—. 

(Hi) (Estimates on p and K = — V~) For I 7 I < A: — 2, 


11021^11^2 <C'eo(l + t)-4 




12V<^|U6 < Ceo(l + t)-4 


-2n-2-M 


and for I 7 I = A: — 1, 

||02i^llL2 <Ceo(l + t)-t-i-l ||02V<^||i6 <Ceo(l + t)-^-'-i 


Remark 1.2 As noted in ]20) . in general, if the shift of the profile is not exactly captured, 
the decay rates for and should be ^ lower than the one given in the above theorem 
even in one space dimensional case. Here, the reason that the above decay estimate holds 
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is that the shift due to the initial perturbation introduced above so that when we apply the 
Green function, the term corresponding to the initial data yields an extra (l + t )“2 decay after 
taking the anti-derivative of the initial perturbation. Moreover, under the condition on 

the initial shift, even though the anti-derivative of the perturbation can not be defined for 
all time as the shift function is not precisely defined, we know that <5* is exactly the final 
shift when t tends to infinity of the profile because the initial perturbation will spread out 
eventually. 

Remark 1.3 We note here that the decay rates of K is higher than that ofV^. However, 
we can only get the decay estimates on derivatives of K up to {k — l)-th order. 


Remark 1.4 Compared with fi,9l 21^, our initial data are the small perturbation of the 
planar waves, instead of the constant states. In the meanwhile, here we can show the 
LI{p E [2, +oo]) convergence rates of the planar waves of the three-dimensional bipolar Euler- 
Poisson equations. This improved the results in n- Moreover, here we only consider the case 
that the far fields of two particles’ velocity in the xi-direction are same, see namely, the 

switch-off case. However, we believe that the same results also hold for the switch-on case. 
Indeed, using the gap function with exponential decay in B we can show the similar results 
for the switch-on case. 


The outline of the proof of the main theorems is as follows. First, we notice that 
the equations for are coupled by V(p, which is expressed by nonlocal Riesz potential 
V<y9 = VA~^K, with K = V~^ — V~. So we need to have some good estimates on K before the 
estimate of V~^ and V~. Luckily we note that K satisfies the damping “Klein-Gordon” type 
equations with an addition good term to perform the energy estimate. The estimates of K and 
V<y9 are given in Section 2, where the algebraic decay rates of K in the L^-norm are derived by 
some delicate energy methods, which will be used to obtain the L'^ijp G [2, -|-oo])-convergence 
rates of the solutions in the subsequent. Next, we use the frequency decomposition based 
energy method introduced in |20] . which combines the approximation Green function and 
energy method, to prove global existence and IP convergence results, see (i) and (ii) in 
Theorem 1.1. This method captures the low frequency component in the approximate Green 
function and avoids the singularity in the high frequency component. That is, we firstly show 
the precise algebraic decay estimate of in the low frequency component, which dominant 
the decay of the perturbation, and then obtain the better decay rates of the high frequency 
component of in the L^-norm by energy methods. For the high frequency component, 
one has an additional Poincare-type inequality to close the energy estimate. Note that the 
lack of Poincare inequality in the whole space is usually the essential difficulty in the energy 
estimate which is in contrast to the problem in a torus. This in some sense illustrates the 
essence of the Green function on the decay rate related to the frequency. 

The rest of the paper is arranged as follows. In Section 2, we will reformulate the system 
around a planar diffusion wave defined in (|1.3I) and then state some known properties of this 
background diffusion wave. In Section 3, we will study the energy estimate of K = — V~ 

and prove part (iii) in Theorem 1.1. The frequency decomposition based energy method will 
be carried out in Section 4 and 5, where in Section 4, we will study the approximate Green 
function and then the main estimates on the low frequency component of and in 
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Section 5, we will study the L? energy estimates on the high frequency component of V^. 
Finally, we will complete the proof to part (i) and (ii) of Theorem 1.1 in Section 6. 


2 Preliminaries. 

In this section, we will first derive the equations for the perturbation functions and 
dehned in (ll.8|) . Then we will recall some results on the background diffusion waves. 


2.1 Reduced system. 

We first derive the system for the perturbation of the nonlinear planar diffusion wave. 
Then, from (11.111 and (|1.4I) . we have the equations for that 

+ (p± + y±)div[/± = V - ([/± • v)(p± + y±) - v^{ufu - hf (2.1) 


where 




P{P^)x,V^ 


+ R^± - Rf ± 


Similarly, the equations for are 

(t/f )i + (p± + v^)-^[Pip^ + F±) - p{p^)U +ut = 4 

p 

and for z = 2,3, 

{ut)t+{p^+v^)-4p{p^+v^)-p{p^)U+ut = 
where 


R-uf = -[uf (xi + 5(a:',t),t)(5z(x',t)]xi, 
Rf= • V(uf + Uf ) + u4{Ut)^,, 
Rf= U^-VUt + ut{Ut)^,,2<i< 


n. 


The equation for p is simply 

Ap = V+ -V- = K, 

it is directly that the perturbed electric field Vp can be expressed by the Riesz potential as 
a nonlocal term 

Vp = VA-^iF. 

Then the system for the perturbation iV^, t/^, p) can be summarized as 

Vt + {P^ + V^)divU^ = Q±, 

iu4)t + + V^)-4V{V^,p^)V^),^ + u4 = Hf± d,^p, 1 < i < 3, 

where V{V^,p^) = [ P'{p^ + eV^)de, and 

Jo 

Q± = - (c/± . V)(p± + F±) - V^iufU - 


( 2 . 2 ) 


(2.3) 


H 


1 - + p±(p±+y±) - > 


h4= - 




P^J^i I ^KP^)^iV- p± r, / / O 

TE -^ p±(p±+y±) -Ri ^ 2 < z < 3. 
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Moreover, we can deduce the equation for V^{x,t) from (12.3^ as 


V,t - A[ViV^,p^)V^] + = QiV^, U^,p^,uf) T div[(p± + y±)V¥p], (2.4) 


where 

g(y±,f/±,p±,h±) 

= [{Rp±)t + Rp±] - (1 + dt){V^ufh, - div[(p± + y±)tC/±] - div[(p± + 
with = {H^, ■ ■ ■ , H^). By linearizing (12.41) around p, we have 
F±-A(a±(x,t)y±) + y± 

= Q{V^,U^,p^,uf) + A{Vi{p,V)V^)TdiY[{p^ + V^)Vip] (2.5) 

=: F±Tdiv[(p± + y±)V</?], 


where a^{x,t) = P'{p^) and 


ri{p^,v^) 



P"{p^ + 92V^)de2)d9i. 


Since 

{Rp±)t + Rp± = {-pf{x,t)6tix',t))^^, 

direct calculation shows that = F{V^, U^,^,u^) is in divergence form, that is. 




( 2 . 6 ) 


where, without confusion, we omit the =h sign, 

F^ = -pt6t - {puiSt)xi, F" =-P{p)xi, ‘2^ < i < n, 
Fii = p{2uiUi + Uf) + V{ui + Uif + Flip, V)V^ 
jpii ^ pii ^ 2[{p + V){ui + Ui)Uii 2<i<n, 

= {p + V)U^Uj + 6,jVi{p, y)y2, 2 < i,j<n. 


Here 6ij is the Kronecker symbol. On the other hand, by linearizing (12.31) 9 around we 
have 


+ {p^)-^V{a^ix,t)V^) + U^ = 


(2.7) 


where, again without confusion, we omit the =h sign, 

H,= Ru - r^iViip, 

H,= -^-^il^^-p-^[F,{pV)V%^-R,, 2<i<n. 


— Ri, 

r2] 
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2.2 Background profile. 


For later use, we include the following known estimates on the background planar wave, 
cf. [28]. By the definition of W{xi), we know for any integer 


2\-N 


sup \ W{xi) - p+\+ sup \ W{xi) - p-\ < C\p+ - P-\{1 + xf) 

a;i>0 xi<0 

|d,\iy(xi)|< C\p+-p.\{l + xl)-^, {h>0). 

Recall that we have assumed in (ll.7p for any multi-index /3, 

First, let us recall the results about the one-dimensional bipolar Euler-Poisson system (11.41) . 

Theorem 2.1 (see w Let {p^,uf){xi,0) be the initial data of one-dimensional bipolar 
Euler-Poisson system (O and fix an integer m >2. If there exists a small positive constant 
Ep such that the initial data {p^,uf){xi,0) satisfy 


\p+- P-\ + 


{p^{z,0) - W{z))dz\\L 2 


+ IIP^('i 0 ) “ ^(■) II + 11^^(■) 0 ) “ V’(') 0)IIh"*+i < Ep, 


and 


/ Xi 

-oo 


(p^iz, 0) - Wiz))dz, + P{WU G L\K) 
then has global classical solution {p^,u^,E) with 

||d'=(p±-iy)(xi,t)||i.(Ri^)<CE,(i + f) 


1 \ fc+i 

T)/ 2 


\\dHu^ + PiW.M^i,t)h,^nu < CEp{l + t) 




1 fc +2 


for any integer k < m -\- 1 if p = 2, and k < m if p = oo. Moreover, there exists a positive 
constant fi such that 

Remark 2.2 Note that {p^,u^) is an intermediate state we constructed to approximate the 
one-dimensional diffusion wave W. The assumptions on the initial data (p^,'ui’^)(xi, 0), with 
p^(xi,0) connecting the two end states p±, can be more regular than the assumptions on the 
initial date of the original problem m- 

Next, from the definition of the planar diffusion waves {p^,vfi) in ()1.6I) . we can readily 

have 

Lemma 2.3 Under the assumptions in Theorem 2.1, the planar diffusion waves {p^,u^) 
defined in satisfy 


sup||<9“(pJ^,u±)(-,x',t)||i2(Ri ) < CEp{l + t)- 

rf! 1 


< CEp{l + t) *'2 


l+|a| 1 

2 4 . 


l+|a| 
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for any multi-index a with |a| < m — 1, and 

Wp'^ - P ll//’"(Ri^) < CEpC 

In addition, for 2 < i <n, 

Wd^^ipti, iui)xi){t)\\L^n-) < CEpe-\ 
ll'9"(p|,(«f)xJ(i)||L-(Rn) < CEpC-K 

Note here again that we can increase the regularity of the assumptions on initial date of 
the one-dimensional problem to get sufficient estimates on the planar diffusion wave. 


3 Estimates on K = V~^ — V . 


In this section, we mainly give the estimate oi K = V~^ — V~ . Recall the linearized 
equation (j2.5p for V^, we see that they are coupled by Vip, which is expressed by the Riesz 
potential as in (12.2p . i.e., Vip = So we need to have some good estimates on K, 

thus V(/j, before the estimate of and V~. To begin with, we give a lemma on the relation 
of V</? and K. 

Lemma 3.1 If K G for any integer I > 1, then V<p G IR^’®(R”). 

Proof. Note that be expressed by the Riesz potential 

V<y9 = = n*K, 

where E = |27r^|“^ thus IZ = ■ Here n = 3 is the space dimension. Then by Hardy- 

Littlewood-Sobolev inequality [24j we have 

||V(^||^6 = ||7^*R||^6<C||R||^2, (3.1) 

and similarly, for any multi index Iqj < k, 

WViphe = \\n*d"^K\\Le<CWK\\L., (3.2) 

that is, V(/? G if RT G ^'(R"). 

Remark 3.2 By Sobolev injection, this lemma automatically indicates 

Vip G R~(R”). 


Now we start to estimate K. The equation for K, from (12.Sp . is 
Ktt - A{a+K) + Kt = E+-E- - div[(p+ + p" + R+ + V-)Vip\ + A[(a+ - a-)R-], 
where E^ on the right hand side are defined in ()2.6p . This equation can also be written as 


Ktt - A(a+R) + Kt + {p+ +p- + V+ + V-)K 
= F+ -E- - V(p+ + p- + V+ + V-) ■ Vip + A[(a+ - a")!/"]. 


(3.3) 
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note here the last term on the left hand side is a good term, which ensures the closure of 
energy estimate for K. Note that the lack of such term in equation of is the main difficulty 
in energy estimate thus we will use the frequency decomposition method introduced in m- 
To proceed, we first give the a priori assumption 


M{t) = 


(3.4) 


max < sup ( 1 +s) 

^ 0<s<t,|a|<fc-2,p>2 


n r-i 1 ^ I |q:| + 1 I 

2 (l-p) + — 


kl + l 


•,s)||lp, sup (l+s)4^ 2 

0<5<t,|Q:|=/c,A:—1 




V^{;s)h2, 


sup (l + s)2' P 

0<s<t,\a\<k-2,p>2 


2 2 


0<s<t,\a\=k,k—l ^ 


Under the assumption in Lemma 2.3 and the above a priori assumption, it is easy to check 
that for any multi-indies a and 7 , the nonlinear terms in satisfy 




< CEpe-^, 


|ct| < k, 


\\d^yF^mLP(Il^)<CM^{l + S, 




(3.5) 


2 ^+ 2 p^ I 7 I < A: - 2. 


Here we should indicate that the above estimates hold for p > 1. 

Now we perform energy estimates. Multiply equation (13.3p by K and integrate the 
resultant equation over R"", we have 


dt 


[ + KKt)dx- [ A{a^K)Kdx+ [ {p+ + p- + + y~)K^dx - [ K^dx 

Jk^ 2 Jjln Jjln J^n 


/ K{F+-F-)dx- V{p^ + p-+V++ V-)-V^KdxF A[(a+- a")U"]iLdx. 

/R'l 7R" 7R" 


Note 


[ -A{a^K)Kdx= [ V(a+iL) • ViLdx = [ o+|ViL 

JR" JR" JR" 

Moreover, using Cauchy-Schwarz’s and Young’s inequality, we have 


dx -- f {Aa^)K'^dx. 

2 Jr" 



v(p+ + p-+V+ + V-)- VifKdx 


< 

< 

< 


\\Vip\\L^\\V{p+ + r + + V-)\\l2\\K\\l2 

C\\V{t + p- + V^ + V-)h2 \\K\\l, 

C{Ep + M{t)){l + t)--^ \\K\\l,, 



a-)V-]Kdx = 

< 

< 



-a-)V-]-VKdx 


C||V[(a+-a-)U-]||i 2 ||ViL||i 2 

CEpe-^^M{t){{l + 1)-?-5 + (1 + t)-f-i)||ViL||i2, 


since the Sobolev norm of a"*" — a has same exponential decay as p'^ — p in Lemma 2.3, 
and further 


/ 

JR" 


KF^dx<eo I K^ + C{eo) I (F^)'^ < eq I + Cieo){E^p + + t) 


p ±\2 


-§n-4 


/R" 


/R" 


'R" 
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here we have used (I3.5p in the above estimate. Combine above estimates and good decay 
properties of in Lemma 2.3, we have 


^ f {lK^ + KtK)dx- [ K^+ [ {K^dx + \VK\^)dx 

dt J^n 2 J^n 


(3.6) 


< C{E^p+M{t)^){l + t)-2 


— Sn—4 


Here and in the subsequent we use the fact that is strictly positive and bounded 

from below. 

Next, multiply equation (|3.3p by Kt and integrate the resultant equality over R”, we 

have 


55' 


/R" 


dx- [ A{a^K)Ktdx+ [ Kfdx+ [ {p+ + p- + V+ + V-)KKt 

JR" JR" JR" 


dx 


/R" 


Kt{F+-F-)dx- V(p+ + p“+ F+ + H“) •V<y9iLtdx+ / A[(a+- a")H"]iLtdx. 


/R" 


/R" 


It is easy to compute 

[ -A{a^K)Ktdx= [ V(a+iL) • ViLjdx 

JR" JR" 

= / a+l-dt\VK\‘^dx+ [ Va+• (iLViLt)dx 

Jr.^ 2 J^n 

= [ a^l-dt\VK\‘^dx + [ ■ dtiKVK)dx - [ ■ {KtVK)dx 

Jr^ 2 J^n J^n 

= -^ [ a'^]-\'VK\‘^dx — [ dta~^]-\VK\'^dx —f ^Aa^K^dx 

dt J^n 2 J^n 2 dt Jjln 2 

+ [ -dtAa^K^dx - [ Va+ • {KtVK)dx, 

JR" 2 jR" 

in which the last term on the right satisfies 

I [ ■ {KtVK)dx\ < C\\Va^\\Loo\\Kt\\L 2 \\VK\\L 2 

JR" 

< ei\\Kt\\L2 + C{ei)Fl{l + t)-l\\VK\\l,. 

Next, note also that p^ + is strictly positive and bounded 

[ (p+ + p- + H+ + V-)KKtdx 

JR" 

= 4/ ]-{p++p-+V++ V-)K^dx- [ ]-dt{p++ p-+ V++ V-)K^dx. 

dt Jun 2 J^n 2 

Also note 

[ v{p+ + r + V+ + v-)-Vip Kt < c\\Vip\\L^\\v{p+ + r + v+ + v-)\\L2\\Kt\\L2 


/R" 


< ^2\\Kt\\\2 + C{£2){Fp + J\A‘^){1 + t) 2||iL||2,2, 


[ A[{a+- a-)V-]Ktdx < e3\\Kt\\l2 + Ciss^F^^ + M^)t-^\ 

JR" 
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and 


/R' 


KtF^dx < £0 [ K^dx + C(so) [ {F^fdx 

J JR" 


< So Kfdx + C{eo)iE^p+M^){l+t )-2 


-§n-4 


/R" 


then combine above estimates to get 
dt 


< CE^ 


[ {Kf + K^ + \VK\^)dx+ [ K^dx 

JR" JR" 

[ {K^+ \VK\^)dx + C{E^p + M^){l + t) 

jR" 


(3.7) 


-|n-4 


Then multiply (j3.6j) by 4 then add to (j3.7p . and if we assnme that Ep + Ai is small enough, 


then 


4 / iK^ + |Vi7p + Kf)dx + [ {K^ + \VK\^ + Kf) < C{eI + 7W^)(1 + t)-i 


-2n-4 


dt 


/R" 


'R" 


by Gronwall’s inequality, we have 



+ \VK\^ + Kf < C{El + M^){1 + 


then 

\\K\\l 2 < C{Ep + M^)il + 

For higher order derivatives (see also the estimate on high frequency part of below), take 
d'^ on both side of equation ()3.3p . multiply by d'^K and d'^Kt and integrate, respectively, 
then combine as above estimates for lower oder derivative. Note that the last term on the 
right hand side of (|3.3p A[(a''' — a~)V~] has already second order derivative and the a priori 
assumption p3.4p has control of derivatives up to k-th order, so we can only carry out the 
computation for derivatives with order I 7 I < k — 2. Similar arguments as above, we have, for 
I 7 I < A: - 2, 

and for I 7 I = A; — 1, 

\\d2KU2<Ceo{l+t)-"^^-^-l (3.9) 

Combine the estimates (j3.1[) - ()3.2p . (I3.8p - (j3.9p and nsing Lemma 3.1, we have the resnlts of 
(iii) in Theorem 1.1. 


4 Approximate Green Function and Estimates on the Low 
Frequency Component. 

In this section, we will give approximation Green function of the equation to as in 
[20] . which is used to get LP estimates on the low frequency component of V^. Recall the 
linearized equation ()2.5p . 

Vtf - A(a±(x,t)F±) + T div[(p± + R±)V^]. (4.1) 
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We slightly abuse notations by dropping ‘=h’ sign without confusion in the following estimates. 
Note the main difference of (j4.ip to the linearized equation in [20] is that we need to consider 
the coupling term ‘=pdiv[(p=*= + V^)Vip\’ in the present setting. 


4.1 Approximate Green Function. 

In this subsection, we study the approximate Green function for (|4.ip . For convenience 
of readers, we briefly repeat the construction of the approximate Green function in below. 

Let G{x, t] y, s) be the approximate Green function for the homogeneous part of (|4.ip 
which meets the basic requirement 


G{x, t] y, t) = 0, Gtix, t] y, t) = 5{x - y), 


where 5 is the Dirac function. Multiplying (14.11) whose variables are now changed to {y, s) by 
G and integrating with respect to y and s over the region R"' x [0, t] to get (note here and 
below we dropped the notation ‘zh’) 

V{x,t)= [ Gs{x,t-,y,0)V{y,0)dy- [ G{x,t]y,0){V+ Vs){y,0)dy 

+ [ [ {Gss - a'^^yG - Gs){x,t;y,s)V{y,s)dyds 

Jo 

n G{x,t;y,s)F{y,s)dyds 

-f f G{x,t-,y,s)dw[{p + V)V(p]{y,s)dyds. 

Jo Jr" 


(4.2) 


+ 


If a{y, s) is a constant and G is the Green function of the homogeneous part of (|4.ip . then the 
third integral in above is zero. However, when a{y, s) is not a constant, it is difficult to give 
an explicit expression of the Green function. Therefore, we will use the approximate Green 
function constructed in [281 120| . The idea is to hrst consider the linear partial differential 
equation 

duV - pAV + Vt = 0, (4.3) 

where /r is a bounded parameter with Gq < fi < Ci, and denote its Green function by 
G^{y;x,t), whose Fourier transform 


G«(/i;e,i) = 


gA+i _ gA_f 

A+-A_ ^ 


where 


A±«) = i(-i±yrvi^). 


Denote & = + E with 

E+ = yoe^+\ E-=r,oe^-\ = {X+- X_)-\ 

The approximate Green function is dehned by 

G(x, t; y, s) = G^{a{y, a{t, s));x -y,t- s), 


(4.4) 
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with a{y,a{t, s)) = P'{p{y,a{t,s))), and the function (T{t,s) is chosen such that cr{t,s) E 
C^{[2, oo] X [0, oo]), 

s, s > t/2, 


a{t,s) = 


t/2, s<t/2-l, 


and 


< C, s e {t/2-l,t/2). 

I<li+l2<3 


Notice that a s) < (7(1 +1) ^ for t > 2 so that we have by Lemma 2.1 
(1 +t)|5^a(y,(T(t,s))| + (1 + t)^|(9^o(?/,cr(Ls))| < CEp, 


(4.5) 


where Ep is defined in Lemma 2.1. 

Notice that the decay of the derivatives of the function a{y, a{t, s)) with respect to time 
will be used in the following analysis. Recall that the approximate Green function defined in 
(14.40 is not symmetric with respect to the variables (x, t) and {y, s). However, straightforward 
calculation gives their relations as 

= -dyfi + da{G^) dtG = -dsG + a,(G«) (a, + at). (4.6) 

Denote the low frequency component in the approximate Green function G{x, t] y, s) by 

Gl{x, t] y, s) = x{Dx)G{x, t; y, s), 

where xi^x), Dx = is the pseudo-differential operator with 

symbol xH) as a smooth cut-off function satisfying 


x(0 = 


1, iei<e, 

0 , |e|> 2 e, 


for some chosen constant e in (0, eo) with £o = ^ ™in |l) upper bound of p 

in (j4.3|] . Moreover, we have 

GLix,t-,y,s) = 72 ^/ xiOe^^'"~'^^^&{a{y,cr{t,s)),C,t - s)d^ 

= G\{a{y,o(t,s))\x-y,t-s). 

Therefore, it is direct to have the following proposition [20j . 

Proposition 4.1 For q € [l,oo] and any indices h,l,a and ( 5 , we have 

sup||(9“(9^(9'af Gl(-, t; y, s)||L9(Rg) <G{l + t- s)“ ^ (i"? 
y 

sup||<9"9^<9i9fGL(x,t;-,s)||i,(R,n) < C{1+ t - 


_ n/-i 1\ 2 + |q:| + |/3| 

2 2 


_ n/ I 1 ^ 2 min{l-\-h,l)-\-\a\ + \0\ 

2 V-*- „ J 2 
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4.2 Estimates on the Low Frequency Component. 

In this subsection, we will establish the estimates on the low frequency component 
by using the approximate Green function. Assume that |a| < A: in this section. To derive the 
U’ estimates for the low frequency part, recall (14.2p . and set 


I? = 


_ 

-'2 — 


T°‘ — 
^3 — 


-'4 — “ 


ra _ 
-'5 — 


x{Dx) [ d^Gs{x,t-,y,0)V{y,0)dy = [ d^{GL)six,t;y,0)V{y,0)dy, 
JR" JR" 

[ d^GLix,t-,y,0){V + Vs){y,0)dy, 

JR" 

n dxRcAx^t-y.s) V{y,s)dyds, 

t," 

n d^Giix, t; y, s)F{y, s)dyds, 

t," 


0 JR" 


5“Gl(x, t; y, s)div[(p + V)Vip\ (y, s)dyds, 


where 

and 

Since 

we have 


Rg = Gss{x,t-,y,s) - a+{y,s)AyG{x,t-,y,s) - Gs{x,t;y,s), 
X{Dx)Rg = RGl- 

(G'i - + G\){a{y,s)\x - - s) = 0, 


RGL{x,t]y,s) 

= G'l; 0 ,o(®( 2 /’ s)]x-y,t- s)asiy, af - ‘^G{.^^^{a{y, s);x-y,t- s)asiy, a) 

+G'l;o(®(2 /> s)-x-y,t- s)assiy, cr) - G{.Q{a{y, s)-x-y,t- s)asiy, a) 

+a{y,s)[ E[G'i-0i(o(y>'5);^ - yR- s)ay,{y,a) 

^i=l 

s)]x-y,t- s){{afy.){y, u)] + G\.^{a{y, s);x-y,t- s)Aya{y, a 
+ [(“(?/>cr) - a(y, s))AG\^{a{y, s)]x-y,t- s)] 

=: Rh,+Rl,. 

Here Rq^, z = 1, 2, is the corresponding term in the above summation in the above equation. 

To denote the derivatives, we use the notations G^.Q(a;x,t) = daG]^{a-,x,t), G^j^.-{a;x,t) = 
dxiG{{a;x,t), G{.^^-^{a;x,t) = dtG{{a; x,t), G{.^{a-,x,t) = dadxiG{{a;x,t), andG{.^.^{a-,x,t) 
dadtG^i{a-, X, t) etc.. 

Then, set XL{x,t) = x{Dx)X{x,t), from above notations we have 
d^XUx, t) = If + If + If + If + If. 

We will estimate the right hand side of above term by term. The terms from If to If are 
similar to the estimates in m- by Proposition 4.1, it is straightforward to obtain 

Wlfhnn^i) < C(l + t)-t(i-i)-^||Eo||Li. 
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For set 

My) = My,^) + ^^(^, 0 ), 

where i' (with ± omitted) is defined in (II.91) . Then 

l-^?l = I [ dyi^xGL{x,t]y,0)My)dy\- 

JR" 

Also by using Propositions 4.1, we have 

77. /-I 1\ |q:| + 1 

l|i^o||Li- 


77 1\ |q:| + 1 


We now turn to estimate the term Jg which is the error coming from the approximate 
Green function. For illustration, we only consider 


= [ [ d^G{.Q{a'^{y,a),x-y,t - s)at{y,a)V{y,s)dyds, 
Jo JR" 


and 



j" = / / d'^RGMq^y^^)^iy^^)dyds, 

to JR" 


because the other terms in q can be estimated similarly. Note that (14.51) gives 

\at{y,(7)\ < CEp{l + tM, 

then we have, for I 7 I < A: — 2 , 

II-A7IIlp(RS) < / 11 / d2G{.Q{a{y,a),x - y,t - s) asiy,a) V{y,s)dy\\LP(RM^ 

Jo JR" 


< CEpM 


ntl2 

/ (1 + t - s) 

JO 


1 ItI 




2 {1 + tM{l + s)~^^^~M^Gds 


+ [ (l+t-s)-t(l-T)(l+t)-l(l + s) 

Jt/2 


I'l—H+i 


2 ^^ 2 dg 


2 nJ 2 


< C.FpA4(l+ t)“2^"-p 

and for I 7 I = A; — 1 and k, 

r-t/2 


II'^7IIlp(R") — CEpAi I (1 +1 — s) 2^ p) 2(1 1 ) 


+CEpM 


E l7l+2-fc 1 1\ fc-2+1 

/ (1 + t — S) 2 (1 -)- t) 1(1 -)- s) 2 ( p) 2 ds 

Jtl2 


nr-i 1\ I 7 I + I 

2 *' D-' 2 


For q, since 


< CEpMil + t)~^'-^~p 


\a{y,s) - a{y,a)\ < // |a^(r/,r)|dr < 


CEpQ{t, s), s < t/2, 
0 , s>t/2, 
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where 


0 (t, s) = (1 + t - s){l + 

and h can be any positive integer. By using Proposition 4.1, we have 

11^2 \\lp{R’^) ^ CEpM. (1 + t —s) 2 p 2 0 (fjs)(i + s) / ds. 
Jo 


By noticing that 

ft/2 


(1 + + sY^I’^ds = (1 + + + s)2-7r 


1/2 




we obtain 


_ILM |q=I + 1 

2 2 


II'^^IIlp{R 5) ^ C'£^pA4(l + t) 2 p 
Thus, combine the above estimate to have 


\m\LPiRt/)<CEpMil + t)-^^^--p^--^. 


Next for recall that F (with ± omitted) satisfies (|3.5p under the a priori assumption 
dsai), then for I 7 I < k — 2 , 


1-^4 lliP(RS) - 


+ 


82Gl Fdy\\LPCRn)ds 

(- E + E F*^)dy||LP(Rg)d5 


/o Jk 
f 2 r 


fO JK 
pi 


/ 2 n /I 1 \ I 7 I + I / - n /-I i \ 171-0-1- 

< CEp (l + t-s)" 2 d-p) —e-^ds + CEp (1 + t - s)" 2 


1 


+CJY^ f \l + t- s) 2 ^^ p) + s)- 0 +i)+tds 

Jo 

+CM‘^ [ (l + t-s)-^(l + s)"^"+^+^^+^ds 


2 ^n , 


< C{Ep + M^){l+t) 

The cases when \'j\ = k — 1 and k can be estimated similarly, and the only difference is the 
estimation on the terms like 

f II [ d2dy^y,GL T*^dy||LP(Rs)d5, l7l < A: - 2. 

4R" 

On the other hand, these terms can be estimated by replacing the derivatives of Gl w.r.t. x 
to the derivatives of w.r.t. y using (14.6p . Then by using integration by parts k — 2 times 
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to transfer the derivatives on Gl to , we have by (4.8) and Proposition 4.1 that 


|7l+2-(fc-2) 


dlGi dy^y^F^^dy\\Lr>in!i)ds < CM^ I (1 +1 - s)"'" ' ■ 2 '' (1 + 5)-(-+i+^)+ids 


/R" 


< C'A42(1+ 


Therefore, we have the U’ estimate on If as 


_'ri(^ 1\ |Q=I + 1 

2 r>) 2 


1-^4 IIlp(R^) ^ (1 + 0 ^ ^ 


For If, we write 


1 " — 
-'5 — 


/o 


+ 



4 JR” 


dfVyGL{x,t;y,s) ■ [{p+ V)VyI\{y,s)dyds 
dfGiix, t; y, s)div[(p + P)V(/?](y, s)dyds 


— • _i_ r" 

— • ^ 5,1 + -'5,2) 


in which the first term satisfies 

t 

Ili<C j" \\dfVyGL\\LA\{p + yWl^\\Fds. 

Jo 

To estimate 

\\{p + V)Vp\\li <G\\p + V\\^,jVip\\L.. 

we note that the Ls norm of p + V{= p, positive) can be controlled by its and norms 
by interpolation, and also note the norm of p (with ± omitted) is conserved because the 
conservation of mass, norm can be bounded by the a priori assumption, then we have 

n 1 \ |q:| + 1 


For the term If 2 : we estimate 

Ii ,2 < C ||5“Gi|Ui||div[(p + P)V^] Uvds. 

Note that both divV^? = K and V(/9 are in L°° thus have good decay properties, then the 
above term decays faster than that of Ifi, then we have 

n /I 1 ^ |q:| + 1 

WI^Wlp <GM^{l + t)~^^^~F 


In summary, by combining all above estimates, we have the estimates on the low 
frequency component of in the following theorem. 

Theorem 4.2. For |a| < k, we have, 

Wd'V^mLp < C{Eo + M^){l + 
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where Eq = max{||yo=^||ii, ||z:^^||z,i, \\{V^,U^)\\Hk, \\V^^{0)\\Hk-i,Ep}. 

As an immediate consequence, we have the estimate on the derivatives of order higher 
than the k-th for the low frequency component because 

\\d^,d^V^{t)\y = = e\\d^V^{t)\\L.. 

Thus, we have the following corollary. 

Corollary 4.2 For any I 7 I > k, we have 

WIIl 2 < C{Eo + M^)e{l + t)-'l^^-h^-^. 


5 Estimates on the High Frequency Component. 


In this section, we will carry out the energy estimates on the high frequency component. 
Recall the linearized equation (12.5h . set xiO = 1 - X(0 and V^{x,t) = xiD^)V^{x,t). By 
taking x{Dx) on both sides of (12.5p and integrating its product with and (V^)t over R"" 
respectively, we have 

4 / V^iV^)tdx-[ {{V^)tfdx- [ Ax{aV^)dx + ^ [ l{V^fdx (5.1) 

ai J-£in J^n dl J A 

= [ V^xF^dx^ [ {V^)txdw[ip^ + V^)Vip], 

JR" JR" 

and 

4 / hvi)tfdx- [ {V^)tAxiaV^)dx+ [ {{V^)t?dx (5.2) 

at Jj^n Z Jj^n J^n 

= [ iV^)txF^dxT [ V^xdiv[{p^ + V^)Vip]. 

Jr" Jr" 

Again, we slightly abuse notations by dropping the ± sign without confusion in the following 
estimates. First for the third term on the left hand side of ()5.ip as follows. That is. 


-[ VHAxiaV)dx=[ a\VVH\‘^dx- [ Vh V[Vx,a]Vd: 
Jr" Jr" Jr" 

where [A, R] = Ao B — B o A denotes the commutator. Since 


(1 + IIVa|h°° + (1 + t)||Aa|h°o CEp, 


where Ep is defined in Theorem 2.1. It is straightforward to show that 



IV[Vx, a]V\‘^dx < CE^pM'^il + ^)■5"^ 


where Ai is defined in (3.4). Thus 


/R" 


Vh V[Vx,a]ydx\ < p 


\VHfdx + CElM\l + t)-^-\ 


/R" 
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We now turn to estimate the second term on the left hand side in (j5.2l) . That is, 


- Lr.{VH)t Ax{aV)dx = / {VVH)t aVx{VH)dx + {VVH)t[Vx,a]Vdx 

[ a\VH\^dx-l [ atlVnl^dx- [ {VH)MVx,a]Vdx, 

JR" ^ JR" JR" 


/R" 

2 dt 


in which, we have 


/ at\V{VH)\^dx\<CE^pil + t)-^ [ |V(Vff)| 2 dx, 
JR" JR" 


and 


/ {VH)MVx,a]Vdx\<rj [ \{VH)tfdx 
JR" JR" 


+ CE^M^{l + t)-2-\ 


For VnxF'dx and J-^riiVH)txPdx on the right hand side of (jS.ip and (I5.2p . by using 
Lemma 2.3 and the definition of Ai in (13.4p . we have 


< CEpe-^, 


Wf^wl^ - 

+ \x\<k-2, 

where F and defined in (12.6p . Further, it is straightforward to check that 

I [ VuxFdxl <7] [ iVnl^dx + Cir]){E^p + 7W^)(e-' + (1 + t)-2(”+i)+t), 

JR" JR" 


and 


\[ {VH)txFdx\ < rj [ \{VH)tfdx + C{rj){E^p+M^){e-^ + {l + t)-^<^^+^F^). 

JR" JR" 

For the term Jj^„(V//)txdiv[(/0 + 'F)V(/?], we have 
[ {yH)tXdiv[{p + V)Vip]dx 

JR" 

: f {VH)tdx + C{e) f {(liY[{p + V)VLp\fdx 

JR" JR" 


< e 

= e 

< e 

< e 


/R" 


{Vnftdx + C{e) / [V(p + V)-\/^+{p + V)Kfdx 


'R" 


j^yH)ldx + 2C{E)(^\V{p + V)\\l4V^\\l^ + \\p^ 

[ {Vnftdx + C{e){El + M^){1 + t)-|^-l“l-2. 

JR" 


The term VHXdxvHp + V)Vip\dx can be estimated similarly. 

To close the energy estimate, one needs the following important fact about the high 
frequency part: 


'R" 


\VVH\^dx > e / \VH\^dx. 


/R" 
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This is a Poincare type inequality which holds only for the high frequency part in the whole 
space. By integrating (j5.ip and (j5.2p over [0, t] and multiplying (15.ip by some suitably chosen 
constant 0 < A < 1, when r] is small, the combination of above estimates give 


< 


L 


{\Vh\^ + \{VH)tf + \VVH\^mdx + f f {\VHf + \{VH)sf + \VVHf)dxds 

[ {\VHf + \{VH)t\^ + \VVH\^mdx + {E^p+M^) [\l + sr^-^ds\, 
jR" Jo 


for some positive Denote 


Ht)= / i\VH\‘^ + \{VH)t\‘^ + \^VH\^)dx. 


/R" 


Then the above inequality gives 

pi 


+ J^{s)ds<C{J^{0) + {El + M^) [ {l + s)-^-^ds). 

Jo Jo 

By using the Gronwall inequality, we have 

J^{t) < Ce-f^\j^{0) + {El + M^) [ e'^"(l + s)"t-3ds). 

Jo 

Hence, we have 

\\VH{tWH. + \\{VH)t{t)\\h 

< e->^^{\\VH{0)\\m + + C{El + M^l + t)-t-3 

< C{El + M^){l + t)-'^-^, 


where Eq is defined in Theorem 4.2. 

Next, we will derive the energy estimates on the higher order derivatives of the high 
frequency component, that is, for 0 < |a| < A: — 1. 

In the rest of this section, we assume 0 < |q;| < /c — 1. The estimation can be obtained by 
induction on lal. Assume that 


jR" 


+ \Vd'^VH\^)dx < C{Ei + M^){1 +1) 


f-(l7l+3) 


holds for any multi-index 7 with I 7 I < |a|, we want to prove 


/R" 


{\d^VH\^ + \d^{VH)t? + \Vd^VH\^)dx < C{El + M^)(l + t)-f-(l“l+3). ( 5 . 3 ) 


Taking d°‘x on (|2.5p . neglecting the ± sign without confusing, and integrating its product 
with d^Vn and d°‘{VH)t over R"" respectively, we have 

4 /’ d^VHd^{VH)tdx- [ \d^{VH)t\^dx- [ d'^Vn Ad^x{aV)dx + ^f hd^Vnl^dx 
at J^n Jjln dt JI 

(5.4) 
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and 


= / d^Vn xd'^Fdx^ [ x5"div[(p± + y^)V(^], 

JR" JR" 

4 / / d'^{VH)t Axd^{aV)dx+ [ \d^{VH)t\^dx (5.5) 

a* JR" ^ JR" JR" 

= / x<9"Fdx^ / 5“(yff)i x5“div[(p± + y±)v^]. 

JR" JR" 

For the third term on the left hand side of (I5.4p . we have 

-/ d^Vn Axd^{aV)dx = ! aid^^VVnfdx - [ d'^Vn V[Vxd'^, a]Vdx. 

JR" JR" JR" 


Since 


it holds that 


IlSfallioc <CEp{l + t)Am^ 




I [ O^Vh V[Vx 5“, aJ'Fdxl < 77 / pdj: + + t)' 

JR" JR" 

Similarly, for the second term on the left hand side of (I5.5h . we have 

-[ d‘^{VH)tAxd^{aV)dx = ^[ ^|V9“VHpdx- / '^\Vd^VH\^dx 

JR" at Jj^n Z Z 

-[ d^{VH)tV[Vd^x,a]Vdx, 

JR" 


where 


/ V[Vxa“, a]ydj;| < r, [ \d^{VH)t\^dx + C^E^pM^l + t) 

JR" JR" 


-f-3-|«| 


For the terms d^^Vnxd^^Fdx and d°‘iyH)txd°‘Fdx on the right hand side of ()5.4n and 
(I53|), we only estimate the second one because the estimation on the first is easier. Notice 
that the estimation on the terms with derivatives of order less or equal to |a| + 1 follows 
directly from the definition of A4 in ()3.4[) . Thus, we consider the terms with derivatives of 
order higher than |q;| + 1. Firstly, by using the expression ()2.5[) for F, we have 

F= Q +/X{Vi{p,V)V^) 

= [{Rp)t + Rp] - (1 + dt){Vui)^, - div((p + V)tU) - div((/i + V)H) + A{Vi{p, V)V‘^). 

(5.6) 

Since ()2.1I) implies 

divl7 = -{p + V)-\Vt + {u + U)-VV+ {U ■ V)p + VdoNu - Rp), (5.7) 

substituting (j5.7|) in (15.6p . by the definition of H, we have 

F = {u + U)- V((u + U) • VV) + A{Vi{p, V)V^) + 77, 
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where TZ denotes the remainder which contains derivatives of U and K with order at most 1. 
Thus, d°^TZ has derivatives with order at most |q;| + 1 (< k). Then 


/ 

JR" 


d^iVHh xd^Fdx 


NI + N 2 + N 3 , 


with 

Ni= [ d^{VH)txd^i{u + U)-V{{u + U)-VV))dx, 

JR" 

N2= [ d^{VH)tXd^A{Vi{p,V)V^)dx, 

JR" 

N3= [ d^{VH)td^ndx. 

J-RP 

For A^i, we have 


Ni = 


[ d^{VH)tHu + U) ■ V((u + U) ■ Vd^V)dx + {•••} 

JR" 

/ d^{VH)t {u + u)- V((u + U) • Vd^VH)dx + f d^{VH)t [(« + U) • Vf,xWVdx + {•••} 

JR" JR" 

f + vd^{VH)fdx - [ {u + u)f vd^{VH){u + U) • vd^{VH)dx 

m JR" ^ JR" 

- [ V(u + [/) • Vd^iVnhiu + U) ■ Vd^{VH)dx + [ d^{VH)t [(« + U) ■ Vf,x]d^Vdx + {•••} 

JR" JR" 

-17 / 7rl(^ + C^)-Va“(yH)Pdx + iVi,i + A^i,2 + iVi,3 + {---}. 
dt Js^n 2 


Here and in the subsequent of this section, we use {• • • } to denote the terms with derivatives 
of order at most |q;| + 1. It is easy to see that 


|iVi,i + iVi,2 + iVi,3 + {• • • }| < C{Eo + M^){1 + t)-f-(l“l+3). 


Then for N 2 , we have 


N 2 = 


d'^iVH)t xd^AiVip,V)V^)dx 


/R" 


/R" 


xVy d^AV + 2d^{VH)t xV{p, V) V d'^AVdx + {•••} 

i 

= : -/V2,l +-^^2,2 + {• • • }• 

By noticing that 

A^ 2 ,i = [ d^{VH)t xV'v d^AV V^dx 

JR" 


'R" 

_± 

dt 


d'^iVH)t {V^V'y) d^AVHdx+ / d^{VH)t [x,V'vy ] d^AVLdx 


/R" 


/R" 


{V'yV^) \Vd^VHVdx + Oi, 
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with 


Oi<T] d‘^{VH)\^dx + C^{Eo + M^){l+t)-'^-^^^\+^'>. 

JR^ 

Similarly, we have 

N2,2 = -2^ [ i'py) |V5“yffpdx + O 2 , 
at J Jin 

with 

02<v [ \d‘^{VH)\^dx + C^{Eo + M^){1 + t)-f-(l“l+3). 

JR" 

We still need to consider the terms from the expansion of 9“div[(p^ + V^)S/(^], 0 < 
\a\ <k — 1 : 

• If all the derivatives cl“div are taken on + V^, it can be bounded by the a priori 
assumption and the fact that Vy? G L°°. 

• If all the derivatives 9“div are taken on V(/9 then 

d^divVip = d^K, 


and recall the good decay properties of d°^K for \a\ < /c — 1 in Section 3, then we have 
better decay on these terms. 


• Other terms can be estimated similarly. 


Again to close the energy estimate, we now use the fact that 


/R" 


\Vxd‘^VH\^dx > e / \d°VH\^dx. 


/R" 


By integrating (j5.4l) and (15.51) over [0, t] and multiplying (j5.4n by some suitably chosen 
constant 0 < A < 1, the combination of above estimates give Therefore, we have 

the following estimates on the high frequency component. 


Theorem 5.1 Under the assumption of Theorem 1.1, we have, for |a| < fc — 1, 

where Eq and M. are defined in Theorem 4-^ and {3-4^ , respeetively. 


6 Proof of Theorem 1.1. 

In the previous two sections, we obtain the following estimates on the low frequency 
component by using the approximate Green function and the high frequency component by 
using the energy method respectively, 

\\d^V^{t)\\Lr <C{Eo + M^){l+t)-^^^--py'-^, \a\<k, ( 6 . 1 ) 

and 

\\d^V^\\Hi + \\d^{V^)t\\L^ <C{Eo + My^){l + t)-^-^, |a|<A:-l. (6.2) 
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It remains to combine (16.ip and ()6.2|) to close the a priori assumption (13.4p . 
Firstly, by taking p = 2 in (|6.ip and combining with (j6.2|) . we have 


\\d^V^{t)\\L2 + + |a| <k. 

Next, by using the Sobolev embedding theorem, from (16.2p . we have, for \a\ < k — 2, 

\\d»v^it)\\L^ < < \\d-v^{t)h2 + ||va“v|(t)||Hi 

Moreover, for n = 3, it holds that —j — < — § — Thus, 

l|5"V^(t)||L- <C{Eo+M^/‘^)il + t)-^-'-^, |al <k-2. 


_ ^ /'I _ 1 \ |q=I+i 

'l 2 r. / 


Then, the interpolation of p6.2p and (|6.3I) leads to 

\\d^V^{t)\\LP < C{Eo + , |a| < k - 2. 

Combining p6.ll) with (16.41) then gives 

||5“F^(t)||LP + + |a| <k-2. 

Now, we turn to estimate by using the equation p2.7p . Note that 

U^{x,t) = e-^U^{x,0) + [ ±V<^){x,s)ds, 

Jo 

and it is easy to check that, for | 7 | < k — 2, 

||a^.ff±(s)||Loo < CiEo + M^)il + s)-^^+^+^\ 

||a^((p±)-iV,(a±y±)(s)|Ucc < C{Ep + Mv){l + 

Wd'TVphoo < Wd'TKWHi < C{Ep + M^){1 + t)-i^ 


(6.3) 


(6.4) 


^-On_2-M 


thus, 


71+2 


(6.5) 


\\dm^it)\\L^ < C{Ep + Mv + M^){l + t)-2-—, |7| <k-2. 

Next, for the L^-norm, we use energy estimate. Multiply (|2.7p by and integrate, 

we have 


L 


p^U^Utdx+ I V{a^V^)-U^dx+ 




/R'* 

Note 


/R" 


7 p^iu^? 


dx = I H^p^U^dxE I p^V(p-U^dx. 

( 6 . 6 ) 


'R" 


/ 

JR" 




T±, 


[ p^U^Ufdx = ^ [ p^iU^ fdx - [ dtp^iU^ fdx 

JR" dt J-£ln J^n 

in which the second term is small, and 

7 V{a^V^)-U^dx\<e [ {U^f + C{s)M‘^{l +1) 

JR" JR" 




(6.7) 
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in which the second term has the decay rate of ||Vy^||^ 2 ) (which is the key observation that 
has better decay than V^). 

The first term on the right hand side of (I6.6p . the nonlinear term, can be estimated as 
in K, which has better decay than (|6.7I) . For the last term with V(p, we have 

\[ p^VcP■U^dx\<\\p^\\L4^^Le\\U^\\L2 < ||p±||i/'||p±||^/'||S|U2||[/±||i2 

JR" 

< e\\U^\\l 2 +C{e)M^ {I+ 

Perform same estimates for 9“(|a| < k), and use the similar argument as for V^, we 
can get 

< M.{1 + , |q;| < A:. (6.8) 

Interpolate (16.51) and (16.81) to have 

I n /-t 1 \ |q:|+2 

WO^U^Wlp + |a| <k-2. 

Then combine the above estimates to get 

Theorem 6.1 Under the assumption of Theorem 1.1, if the initial data (y^,U^) satisfies 
that 


\P+ — P-\ + ll^^(■)0)llL2nLi + l|z^t^(')0 )IIl2 + II^o'^IIr»=uli + ||F)'^(0)||//fc-i + IlF'if 11//*: < eo 


±1 


r±l 


r±i 


where eo > 0 is a small constant, then there exists a unique global classical solution {V^, U^) E 
(^([O,oo),n (^^((O,oo ),to ()2.3p . Moreover, we have 


^2V^\\lp< C(l + rt(^-7)-^, H<k-2, 


\\d2V^\\L2< C{l + t) 4 |7| = A:-1, A:, 


WdfU^Lp^ C(l + t)"5(^-7)-^, |r|<A:-2 


\\d'lU^\\L2< C{l + t) 4 \-f\=k-l,k. 


This theorem implies (j3.4p then closed the a priori assumption, and then it yields the main 
results (i) and (ii) in Theorem 1.1. 

Before concluding this paper, we point out that even though the above discussion is for 
the space dimension n = 3, other higher dimensional cases can be considered similarly. 
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